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Abstract 

Uplink synchronization in orthogonal frequency-division multiple-access (OFDMA) systems is a 
challenging task. In IEEE 802.16-based networks, users that intend to establish a communication link 
with the base station must go through a synchronization procedure called Initial Ranging (IR). Existing 
IR schemes aim at estimating the timing offsets and power levels of ranging subscriber stations (RSSs) 
without considering possible frequency misalignments between the received uplink signals and the 
base station local reference. In this work, we present a novel IR scheme for OFDMA systems where 
carrier frequency offsets, timing errors and power levels are estimated for all RSSs in a decoupled 
fashion. The proposed frequency estimator is based on a subspace decomposition approach, while timing 
recovery is accomplished by measuring the phase shift between the users'channel responses over adjacent 
subcarriers. Computer simulations are employed to assess the effectiveness of the proposed solution and 
to make comparisons with existing alternatives. 
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I. Introduction 

The demand for high data rates in wireless communications has led to a strong interest in 
multicarrier modulation techniques, and particularly in orthogonal frequency-division multiple- 
access (OFDM A), which has become part of the IEEE 802.16 family of standards for wireless 
metropolitan area networks (WMANs) [1]. 

Despite its many appealing features, OFDMA is extremely sensitive to timing errors and carrier 
frequency offsets (CFOs). The former give rise to interblock interference (IBI), while the latter 
produce interchannel interference (ICI) as well as multiple access interference (MAI). To cope 
with such impairments, the IEEE 802.16 standards specify a synchronization procedure called 
Initial Ranging (IR) by which users adjust their transmission parameters so that uplink signals 
arrive at the base station (BS) synchronously and with approximately the same power level. In 
its basic form, the IR process develops through the following steps. First of all, each ranging 
subscriber station (RSS) computes frequency and timing estimates on the basis of a downlink 
control channel. The estimated parameters are used in the subsequent uplink phase, during which 
each RSS transmits a randomly chosen code over a ranging time-slot. As a consequence of the 
different users' positions within the cell, uplink signals arrive at the BS at different time instants. 
Furthermore, since the ranging code is randomly selected, several users may collide over a same 
time-slot. After identifying colliding codes and extracting timing and power information, the 
BS will broadcast a response message indicating which codes have been detected and giving 
instructions for timing and power adjustment. 

From the above discussion, the main functions of the BS during the ranging process may be 
classified as multiuser code detection and multiuser timing/power estimation. Some methods to 
accomplish these tasks were originally suggested in [2] and [3]. In these works, a long pseudo- 
noise (PN) code is transmitted by each RSS over all available ranging subcarriers. Code detection 
and timing recovery is then accomplished on the basis of suitable correlations computed in either 
the frequency or time domains. This approach requires huge computational complexity since one 
correlation must be evaluated for each possible ranging code and hypothesized timing offset. 
Moreover, in the presence of multipath distortions ranging subcarriers are subject to different 
attenuations and phase shifts, thereby leading to a loss of the code orthogonality. This gives rise 
to MAI, which severely degrades the system performance. Alternative solutions can be found 



March 17, 2009 



DRAFT 



2 



in [4] and [5]. In particular, the method in [4] replaces the PN ranging codes with a set of 
modified generalized chirp-like (GCL) sequences and mitigates the effects of channel distortion 
through differential detection of the ranging signals. Unfortunately, this approach is still plagued 
by significant MAI. The scheme discussed in [5] aims at reducing the system complexity by 
breaking the multiparameter estimation problem into a number of successive steps. However, it 
is specifically devised for flat channels and fails in the presence of multipath distortions. 

All previously discussed methods exhibit poor performance in the presence of frequency 
selective channels. Such drawback is partly mitigated in [6] by employing ranging subchannels 
composed by a small set of adjacent subcarriers over which the channel gains do not vary 
significantly. In order to achieve multiuser and multiantenna diversity gains, each RSS exploits 
channel estimates obtained during the downlink phase to select the best subchannel (i.e., the one 
characterized by the largest power gain). Code detection is then accomplished by correlating 
the received frequency-domain samples with the corresponding code sequence and comparing 
the result with a pre-assigned threshold. In the presence of timing offsets, however, the received 
codes are affected by different linear phase shifts and loose their orthogonality, thereby leading to 
residual MAI. A signal design which is robust to multipath distortions is proposed in [7], where 
ranging signals are divided into several groups and each group is transmitted over an exclusive 
set of subcarriers with a specific timing delay. This approach leads to a significant reduction of 
MAI as signals of different groups are perfectly separable in either the frequency or time domain. 
Better results are obtained with the orthogonal signal design presented in [8]. In this scheme, each 
RSS selects a ranging subchannel composed of a specified number of subcarriers and transmits a 
randomly chosen code over adjacent OFDM A symbols. Spreading is thus performed in the time 
domain as the same code is transmitted in parallel over all selected subcarriers. In a perfectly 
frequency synchronized scenario, codes transmitted on different subcarriers remain disjoint at the 
receiver. Furthermore, if the channel response keeps constant during the overall ranging period, 
codes received over the same subcarrier are still orthogonal and can easily be separated at the 
BS. After multiuser code detection, timing information is eventually acquired in [8] through an 
iterative procedure, which exploits the autocorrelation properties induced by the cyclic prefix 
(CP). In spite of the improved robustness against channel distortions, spreading across adjacent 
symbols increases the sensitivity of the system to residual CFOs. In [8] it is assumed that 
during the ranging process frequency errors are so small that the demodulated signals incur only 
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negligible phase rotations over one OFDMA symbol. However, phase rotations may become 
significant if the ranging period spans several adjacent symbols. In such a case, the received 
ranging signals are no longer orthogonal and CFO compensation is necessary to avoid a serious 
degradation of the system performance. 

In the present work we adopt the orthogonal signal design of [8] and propose a novel 
ranging method for OFDMA networks that is robust to frequency errors. To avoid complex 
multidimensional optimizations, we adopt a multistage approach where the number of active 
codes is first determined by resorting to the minimum description length (MDL) principle [9], and 
the multiple signal classification (MUSIC) [10] algorithm is next employed for code identification 
and CFO estimation. Frequency estimates are then used in the third step, where timing and power 
level estimation is accomplished in an ad-hoc fashion. 

The rest of the paper is organized as follows. Section II describes the adopted signal model. 
In Section III we address the problem of code identification and CFO recovery, while Section 
IV is devoted to timing and power estimation. Section V illustrates a method to detect possible 
collisions between RSSs that use the same code and ranging subchannel. Simulation results are 
presented in Section VI and some conclusions are drawn in Section Vn. 

Notation: Matrices and vectors are denoted by boldface letters, with I^v and Oat being the 
identity and null matrices of order N, respectively. A = diag{a(n) ; n — 1,2, ... , N} denotes 
an N X N diagonal matrix with entries a{n) along its main diagonal, while and tr{B} are 
the inverse and the trace of a square matrix B. We use E{-}, (■)*, (■)^ and (■)^ for expectation, 
complex conjugation, transposition and Hermitian transposition, respectively. The notation ||-|| 
represents the Euclidean norm of the enclosed vector while | • | stands for the modulus. Finally, 
round{x} indicates the integer closest to x and [■]k,e denotes the (A;,^)th entry of the enclosed 
matrix. 

II. System description and signal model 

A. System description 

The investigated OFDMA system employs N subcarriers with frequency spacing A/ and 
indices in the set J" = {0, 1, . . . , A?^ — 1}. To avoid aliasing problem at the receiver, a number A^o 
of null subcarriers are placed at both edges of the signal spectrum. We denote by R the number of 
ranging subchannels and assume that each of them is divided into Q subbands uniformly spaced 
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over the signal bandwidth at a distance {Nu/Q)Af from each other, where Nu — N — 2Nq is 
the number of modulated subcarriers. A given subband is called a tile and is composed by a set 
of V adjacent subcarriers. The subcarrier indices in the gth tile (g = 0, 1, . . . , Q — 1) of the rth 
subchannel {r — 0,1, . . . , R — 1) are collected into a set Jq^^ — {i^q]y}]^=o with entries 



The rth subchannel is thus composed of subcarriers with indices taken from J^'^^ = U^Jq J?' , 
while a total of Nji = QVR ranging subcarriers are available with indices in the set JT^ = 
[J^~^ J'^''^'\ The remaining Nd = Nu — Nr subcarriers are used for data transmission and are 
assigned to data subscriber stations (DSSs) which have successfully completed their IR process 
at an earlier stage. We denote by M the number of OFDM symbols in a ranging time-slot. In 
the sequel, we assume that M is a power of two. 

The proposed ranging process develops through the following steps: 

1) Each RSS selects one of the R available ranging subchannels according to some specified 
criterion. In low mobility applications, channel estimates obtained during the downlink slot can 
be exploited to choose the less attenuated subchannel [6]. Here, we follow a simpler approach 
in which the ranging subchannel is randomly selected by each RSS. 

2) The RSS transmits a randomly chosen code of length M during the ranging time-slot. 
Similarly to [8], such a code is transmitted in parallel over all subcarriers belonging to the 
selected subchannel and is taken from an orthogonal set C = {ci, C2, . . . , cm} (e.g., a Walsh- 
Hadamard set) with = [ck{0), Cfe(l), . . . , Cfc(M — 1)]-^ and |cfc(m)| = 1 for < m < Af — 1. 

3) On the basis of the received uplink signals, the BS determines which codes are actually 
being employed and extracts the corresponding frequency, timing and power information. It also 
detects possible collisions between RSSs that use the same code and ranging subchannel. 

4) Once the above operations have been successfuly completed, the BS will broadcast a 
response message by which the detected RSSs can adjust their synchronization parameters. 
Since users that do not find their ranging information in the response message must re-initiate 
the ranging procedure in the next frame, there is no need to notify the collision status to collided 
RSSs. 

Without loss of generality, in the ensuing discussion we concentrate on the rth subchannel 
and assume that it has been selected by K^'^^ RSSs. To simplify the notation, the subchannel 
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index is dropped henceforth. 

The waveform transmitted by the kth RSS propagates through a multipath channel with 
impulse response = [hk{0) , hk{l) hk{Lk — 1)]^. As the channel order Lk is usually 
unknown, in practice we replace by an L-dimensional vector h'^ = [h^, 0, . . . , 0]^, where 
L > max{Lfc} is a design parameter that depends on the maximum expected channel delay 

k 

spread. At the BS, the received samples are not synchronized with the local references. We 
denote by 9^ the timing error of the A:th RSS expressed in sampling intervals while denotes 
the frequency offset normalized by the subcarrier spacing. As explained in [11], users that 
intend to access the network compute initial frequency and timing estimates on the basis of a 
downlink control signal broadcast by the BS. The estimated parameters are then employed by 
each RSS as synchronization references for the uplink ranging transmission. This means that 
during IR the CFOs are only induced by Doppler shifts and/or downlink estimation errors and, 
in consequence, they will be quite small. As an example, consider the IEEE 802.16 standard 
for WMANs with subcarrier spacing A/ = 11.16 kHz and carrier frequency /o = 2.5 GHz. 
In case of ideal downlink synchronization, the maximum CFO in the uplink is 2foVk/c, where 
Vk denotes the terminal speed while c = 3 x 10^ m/s is the speed of light. Letting Vk — 120 
km/h yields Ek < 0.05, which means that the normalized CFO lies within 5% of the subcarrier 
spacing. Timing errors depend on the distances of the RSSs from the BS and their maximum 
value corresponds to the round trip propagation delay for a user located at the cell boundary 
[11]. This parameter is known and given by O^g,^ = 2Rc/{cTs), where is the cell radius 
and Ts = l/(iVA/) the sampling period. A simple way to counteract the effects of timing 
errors relies on the use of sufficiently long CPs comprising Nq > ^max + L sampling intervals. 
This leads to a quasi-synchronous network in which timing errors do not produce any IBI and 
only appear as phase shifts at the output of the receive discrete Fourier transform (DFT) unit. 
Although such a solution is normally adopted during IR, the CP of data symbols should be made 
just greater than the channel length to minimize unnecessary overhead. It follows that accurate 
timing estimates must be obtained during the ranging period in order to avoid IBI over the data 
section of the frame. 



March 17, 2009 



DRAFT 



6 



B. Signal model 

We denote by Y^{iq^i,) the DFT output over the iq^Jh subcarrier of the mth OFDM symbol. 
Assuming for simpUcity that the DSSs are perfectly aligned to the BS references, their signals 
do not contribute to Ym{iq,v)- In contrast, the presence of uncompensated CFOs destroys or- 
thogonality among ranging signals and gives rise to ICI. On the other hand, recalling that the 
CFOs are confined within a small fraction of the subcarrier spacing, the demodulated signals 
incur negligible phase rotations over one OFDM symbol and the resulting ICI can reasonably 
be neglected. Under the above assumptions, we may write 

K 
k=l 

where Nt — N + Nq denotes the duration of the cyclically extended OFDM symbols and 
T^m{iq,v) accounts for Gaussian noise with zero mean and variance cr^. In addition, we have 
defined 

where 

L-l 

Huii,,.) ^Y.^k{^)e-^^^''^-"' (4) 

is the k\h channel frequency response over the iq,v^ subcarrier. The power received from the 
k\h RSS over the ranging subcarriers is defined as 

Q-lV-l 
^ q=Q v=0 

In the following sections we show how the quantities {Ym{iq,u)} can be exploited to identify 
the active codes and estimate the corresponding CFOs, timing errors and power levels. As 
anticipated, in doing so we adopt a multistage procedure in which frequency estimates are 
preliminarily obtained and are next used for timing and power estimation purposes. For the time 
being, we let X < M— 1 and assume that the active RSSs use different codes over the considered 
subchannel [6], [8]. A method to identify possible collisions between RSSs employing the same 
code is presented in Sect. V. 
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III. Code detection and CFO estimation 

We define an M— dimensional vector Y{iq^i^) = [YQ{iq^^),Yi{iq^,y), . . . ,YM^i{iq,,y)]'^ which 
collects the iq^„th DFT output across the ranging time-slot. Then, from ^ we have 

K 

= ^ Sk{9k, iq^u)T{£k)ck + n{iq^^) (6) 

k=l 

where 

T{ek) = cliag{e^'2™"'=^^/^ ; m = 0, 1, . . . , M - 1} (7) 

while n{iq^iy) = [no{iq^^),ni{iq^iy), . . . , nM-i{iq,u)]'^ is Gaussian distributed with zero-mean and 
covariance matrix ct'^Im- Letting e= [ei, €2, ■ ■ ■ , Sk]^ and 6 = [9i,92, ■ ■ ■ , 9k]^, we may rewrite 
Y{iq^,y) in a more compact form as 

Y{tq^,) = C{e)Sie,tq,,) +n{iq^,) (8) 
where C(£) is the following M x K matrix 

C(£) = [r(£i)ci T{e2)c2 ■ ■ ■ r{eK)cK] (9) 

and S{0,iq^^) = [Si{9i,iq^^), S2{92,iq,i.), ■ ■ ■ , SK{OK,iq,u)f ■ From Q we see that Y(zg^^) is a 
superposition of frequency-rotated codes {T(ek)ck}f=i embedded in white Gaussian noise and 
with random amplitudes {Sk{9k,iq,u)}k=i- This model has the same structure as measurements 
of multiple uncorrected sources from an array of sensors. We can thus identify the active codes 
and their corresponding CFOs by applying subspace-based methods. 

A. Determination of the number of active codes 

The first problem is to determine the number of received codes over the considered ranging 
subchannel. This task can be accomplished through the eigenvalue decomposition (EVD) of the 
correlation matrix Ry =E{Y(ig j,)Y^(ig j,)}. In practice, however, Ry is not available at the 
receiver. One common approach is to use the sample correlation matrix 

^ V-lQ-l 

^^ = m7EEY(MY^(^...) (10) 

^ y=0 5=0 
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which provides an unbiased and consistent estimate of Ry. Performing the EVD on Ry and 
arranging the corresponding eigenvalues Ai > A2 > • • ■ > \m in non-increasing order, we can 
find an estimate K of the number of active codes through information-theoretic criteria. Two 
prominent solutions in this sense are based on the Akaike and MDL criteria. Here, we adopt the 
MDL approach which looks for the minimum of the following objective function [9] 

J^{K) = ^k{2M - k) In(Qy) - QV{M - k) In p{k) (11) 

with 




p{K) = . (12) 

i=K+l 

Extensive numerical simulations indicate that a better estimate of K is obtained by replacing 
the smallest eigenvalue Xm with an estimate of the noise power. The latter can be obtained 
in various ways. One possible approach is based on the use of null subcarriers placed at the 
spectrum edges and reads 

M-l No-l 



where Vm(n) is the DFT output corresponding to the nth subcarrier of the mth OFDMA symbol. 



B. CFO estimation and code detection 

Inspection of ^ reveals that the observation space can be decomposed into a signal subspace 
Sg spanned by the rotated codes {T{ek)ck]^^i plus a noise subspace 5„. Since 5„ is the 
orthogonal complement of Ss, each vector in Ss is orthogonal to any other vector in 5„. For 
simplicity, we assume that the number of active codes has been correctly estimated and denote by 
{ui, U2, . . . , ua/} the eigenvectors of Ry corresponding to the ordered eigenvalues Ai > A2 > 
• ■ ■ > Xm- The MUSIC algorithm relies on the fact that the eigenvectors {u/c+i, ue-+2, • • • , um} 
associated with the M — K smallest eigenvalues of Ry form an estimated basis of iS„ and, 
accordingly, they are approximately orthogonal to all vectors in the signal space [10]. An estimate 



March 17, 2009 



DRAFT 



9 



of 6k is thus obtained by minimizing the projection of r{e)ck onto the subspace spanned by the 
columns of = [u^+i ■ ■ ■ ^m]- This leads to the following estimation algorithm 

4 = argmax{^fc(e')} (14) 

e 

with 



^k{e) = 2- (15) 

u^r(5)c. 

In principle, CFO recovery must only be accomplished for the active codes. However, since at 
this stage the BS has no knowledge as to which codes are being employed over the considered 
subchannel, frequency estimates {ii,i2, ■ ■ ■ ,iAi} must be evaluated for the complete set C. 
Next, the problem arises of how to identify the codes that are actually active. The proposed 
identification algorithm looks for the K largest values in the set {'^e{i£)}fL-^^, say {'^u^,iSuk)}k=v 
and declares as active the corresponding codes {c„^}|L]^. The CFO estimates are eventually found 
as Su = [iui,£u2y ■ ■ ■ ,^uk] ■ III sequel, we refer to (1141) as the MUSIC-based frequency 
estimator (MFE), while the described identification algorithm is called the MUSIC-based code 
detector (MCD). 



C. Remarks 

1) A fundamental assumption behind the MUSIC estimator is that the dimension of iS„ is at 
least unitary. This implies K < M, which explains why the number of RSSs over the same 
ranging subchannel should not exceed M — 1 . 

2) Let K < M — 1 and assume that two or more RSSs share the same ranging code over the 
considered subchannel. In such a case, although the MDL can still provide the exact number 
of active RSSs, the MCD is not capable of identifying the corresponding codes as it implicitly 
assumes that any given code is used by no more than one RSS. As an example, let M = 4 
and K = 3 with one RSS employing code Ci and the other two RSSs C2. In this situation, it is 
likely that the MCD declares as active Ci and C2 plus a third code (C3 or C4) which is actually 
turned off and corresponds to a phantom RSS. In a similar way, the MFE will provide three CFO 
estimates, two of which are associated with active users while the remaining one corresponds 
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to the phantom RSS. It is worth observing that, when K = K, the presence of a phantom RSS 
always implies the mis-detection of an active user, which is referred to as undetected RSS. 

3) In multiple frequency estimation problems the ESPRIT [12] represents a valid alternative 
to the MUSIC [13]. The main advantage of ESPRIT is that it provides estimates of the signal 
parameters in closed-form without requiring any time consuming grid-search. A basic assumption 
behind this technique is the rotational invariance property of the observation vectors, which is 
guaranteed in the presence of complex exponentials in noise. Unfortunately, in general, the rotated 
codes {T{ek)Ck}^^^ in general do not satisfy the invariance property. In this case, the ESPRIT 
cannot be used and we must apply the MUSIC or alternative reduced-complexity schemes like 
root-MUSIC or the minimum-norm method [14]. 

4) A key issue is the maximum CEO that the MEE can handle. To simplify the analysis, 
assume that the ranging codes belong to the Eourier basis of order M, i.e., 

c;.(m) = e^'2™(*^-i)/^, < m < M - 1 (16) 

for A; = 1,2,..., M. Then, in Appendix it is shown that identifiability of the rotated codes 
{T{ek)Ck]^^i is guaranteed as long as the normalized CEOs are smaller than N/{2MNt) in 
magnitude. On the other hand, assume that Sk = N/ (2MNt) and Ek+i = —N/ (2MNt). In such a 
case, the received codes T(ek)ck and T{ek+i)ck+i are identical and there is no way to distinguish 
among them. The above facts say that the acquisition range of MEE is \£k\ < N / {2MNt)- 
Although this result cannot be easily extended to other code designs, extensive simulations 
indicate that it also applies to Walsh-Hadamard binary codes. 

5) The accuracy of MEE can be assessed with the methods developed in [15]. Specifically, 
at high SNR values and for large data records (i.e., large values of QV), the CEO estimation 
errors are nearly Gaussian distributed with zero mean and variances 

- ^'^n - STT^gViVlP, ■ df(e,)C^(s)d,(e,) ^^^^ 

where dk{ek) is an M— dimensional vector with entries {mcfc(m)e-'^'^'"'^'=^^/^}^ro ^'^'^ ^'^{^) = 
Im-C(s)[C^(£)C(£)]-iC^(6). 

6) It is interesting to assess the computational requirement of MEE and MCD. Computing 
Ry in (flOl) approximately involves 2QVM'^ real multiplications plus the same number of real 
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additions. In writing these figures we have borne in mind that Ry is Hermitian and, accordingly, 
only the entries [Ry]m^ ,„2 with mi > m2 need be computed. Performing the EVD on a real- 
valued M X M matrix needs 12M^ floating operations (flops) [16]. For the complex-valued 
matrix Ry, we approximate the flop count as 72 by pessimistically treating every operation 
as it were a complex multiplication. Finally, evaluating the metric "^ki^) in (fT5l) requires (4M + 
2)(M — K) real products plus AM{M — K) real additions for each e. Denoting by the 
number of candidate values e and observing that the maximization in (fT4l) must be performed 
over the entire code set {ci, C2, . . . , Cm}, it follows that a total of 2N^M{AM + 1)(M - K)) 
operations are required to find £„. The overall number of flops required by MFE and MCD in 
the considered ranging subchannel are summarized in the first row of Tab. I. 

IV. Estimation of the timing delays and power levels 

After code detection and CFO recovery, the BS must acquire information about the timing 
delays and power levels of all ranging signals. This problem is addressed by resorting to a two- 
step procedure in which estimates of S(0,Zq j,) are firstly computed over all ranging subcarriers 
and are subsequently exploited to obtain the timing errors and power levels. To simplify the 
notation, in the following derivations the indices {uk}k=i of the detected codes are relabeled 
according to the map ^ k. 

A. Timing recovery 

Let K = K and assume that the accuracy of the CFO estimates is such that e ^ £. Then, the 
maximum likelihood (ML) estimate of S{6,iq^^) based on the model © is found to be 

S(z,,,) = [C^(£)C(£)]-iC^(£)Y(z,,,) (18) 

where C(s) = [r(ei)ci T (62)02 ■ ■ ■ r(e^)c^] and we have omitted the functional dependence 
of S{0,iq^^) on 0. Substituting ^ into ^ yields 

S(«,,,) = S(2,,,) + V(^g,.) (19) 

where ilyiiq,^) = [i>i{iq,u)-,'4'2{'i'q,u)-, ■ ■ ■ li'Kiiq^u)]^ ^ zero-mean disturbance vector with co- 
variance matrix cr^[C^(£)C(£)]~^. From (fT9l) and ([3]) it follows that the entries of S{iq^y) can 
be modeled as 
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Sk{tq,u) = e"^^^'>''"'-"'HkK,) + l<k<K (20) 

for any g G {0, 1, . . . , Q — 1} and u E {0,1, . . . ,V — 1}. The above equation indicates that, 
under the assumptions of ideal CFO recovery and code detection, Sk{iq,v) is only contributed 
by the fcth RSS. Hence, we can use the quantities {Sk{iq,u)} to get timing estimates individually 
for each RSS. 

To see how this comes about, we recall that in practical multicarrier systems the channel gains 
over contiguous subcarriers are highly correlated and their values are almost the same. Hence, 
letting Hk{iq^y — 1) ~ Hk{iq^^) and neglecting for simplicity the noise contribution, from (|20l) 
we have 



The above result indicates that a timing estimate can be obtained in the form 



(21) 



9h = round i — arg 

2n 



■Q-iy-i 

'^^Sk{iq,u - l)^fc(V) 

q=Q u=l 



(22) 



where the round operation is justified by the fact that 6'^. is integer- valued. 

As discussed in [11], IBI is present during the data section of the frame whenever the timing 
error 9k — 0^ lies outside the interval J^e = [L — Nqd — 1? 0], where Ng,d is the CP length of the 
OFDMA symbols during the data section of the frame. Intuitively speaking, the probability of 
occurrence of IBI can be reduced by shifting the expected value of the timing error toward the 
middle point of Ja6», which is given by (L — Nqd — 1)/2. This leads to a refined timing estimate 
91 = 9k — fik + {L — Ngd — l)/2, with fik =E{6'fc — 9k}. Extensive simulations indicate that 
Hk equals the mean delay associated to h^. This parameter is generally unknown at the receiver, 
but can roughly be approximated by (L — l)/2. Combining the above results with ((22)) produces 
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if) 



round <( ^ arg 



■Q-l V-l 

.g=0 U=l 



N, 



GD 



(23) 



which is referred to as the ad-hoc timing estimator (AHTE). 

An alternative approach to get timing estimates from the quantities {Sk{iq,u)} is illustrated 
in [17] by resorting to the least-squares criterion. Compared to this solution, AHTE is simpler 
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to implement and can handle a maximum timing offset of N/2, whereas in [17] the estimation 
range depends on Q and L and is typically much smaller than N/2. 

B. Power level estimation 

Collecting (O and (|20l) . we observe that Sk{iq,u) is an unbiased estimate of Sk{iq,i) with 
variance 

al = a''-[T-\e)]k,u (24) 

where T(6) = C^(£)C(£:). An unbiased estimate of \Sk{iq,u)f' is thus given by \Sk{iq,u)\'^ ~(^k- 
This fact, together with ([5]) and (|24)) . leads to the following ad-hoc power estimator (AHPE) 

Q-iy-i 

^ q=0 u=0 

where is an estimate of a^. 

In case of ideal frequency and noise power estimation (i.e., e = e and a'^ = a^), it can be 
shown that AHPE provides unbiased estimates with variance 

varm} = ^^^^^^ for A: = 1,2,...,K. (26) 

C. Complexity issues 

In assessing the computational load of AHTE and AHPE, it is useful to distinguish between 
a first stage, leading to vectors S(ig j^) in (fTSi) , and a second stage where these vectors are used 
to compute the timing and power estimates. We begin by observing that S{iq^u) is the solution 
of a linear system {£)C{£)S{iq^i^) = {e)Y{ig^i,). Assuming that the entries of C{s) are 
available, computing the Hermitian matrix {£)C(e) involves AK^M flops while evaluating 
{e)Y{ig^iy) needs 8KM flops for any pair (g, z/). Each linear system is efficiently solved 
with approximately 2K^{K + 6) flops by exploiting the Cholesky factorization of C^(£')C(£) 
[16]. A total of AK^M + 2KQV{K^ + 6K + AM) flops are thus required to obtain S{iq,^) for 
g G {0, 1, . . . , Q — 1} and u E {0,1, . . . ,V — 1}. In the second stage of the estimation process, 
timing and power estimates are obtained through (l23l) and (l25l) . respectively. Computing 9k from 
(1231) involves 8Q{V — 1) flops for any active RSS, while AQV additional flops are needed to 
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evaluate in (l25l) . The overall requirement of AHTE and APHE in the considered ranging 
subchannel is summarized in the second row of Tab. 1. In writing these figures we have not 
considered the noise power estimate (fT3l) . which is computed off-line with negligible complexity. 

V. Collision detection 

So far, we have neglected possible collisions between RSSs that choose the same ranging 
opportunity. This implies that each subchannel is accessed by no more than M — 1 RSSs which 
employ different codes. Although a proper system design can reduce the risk of a collision, such 
event occurs with some non-zero probability. As mentioned previously, if i^' > M the MFE 
cannot work properly since in this case the noise subspace reduces to the null vector. Another 
critical situation is represented by the presence of pairs of phantom and undetected RSSs, which 
is likely to occur when K < M — 1 and the same code is shared by two or more RSSs. If the 
undetected RSS employs the same code of a detected RSS, the former will adjust its transmission 
parameters according to the response message transmitted to the latter. Such adjustment may 
have detrimental effects as it is based on incorrect synch information. On the other hand, the 
presence of a phantom RSS is not a big problem as the code associated to the corresponding 
response message will not recognized by any of the active RSSs. The above discussion indicates 
that, although it is not strictly necessary to notify the collision status to collided RSSs, collision 
events must be detected to avoid that incorrect synch information be trasmitted to the active 
RSSs. For this purpose, we consider the following vectors 

AY(«,,,) = Y(«,,,) - C(6)S(z,,,) (27) 

for g G {0, 1, . . . , Q — 1} and z/ G {0,1, . . . ,V — 1}. From ([8]), we conjecture that the entries of 
AY(ig J,) will have small amplitudes if C(£) and S{iq^u) are good approximations of C(£) and 
S{iq,u), respectively. Such situation occurs when the estimates K and e are sufficiently accurate. 
Indeed, using standard computations it can be shown that 

E{|| AY(^,,,)||'} = a\M-k) + 6{e, e, z,,,) (28) 

where 
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S{e, e, i 



if fsT = and £ = £ 

2 (29) 
otherwise 



Z^(£)C(£)S(Z,,,) 

and Z(£)Z^(£) is the Cholesky factorization of the MxM matrix Im-C{,£)[C^ {,£)C{£)]-^C" {e) 
Inspection of (|29l) reveals that 5{£,£^iq^y) is minimum when and £ are perfectly estimated, 
while much larger values are expected in the presence of colliding RSSs due to the poor quality 
of K and £. This suggests the use of 5{£,£,iq^y) as a collision detection metric. Unfortunately, 
this quantity is unknown and must be estimated in some manner. Replacing the expectation of 
||AY(im,)||^ in (l28l) with the corresponding ensemble average, an estimate of 5{£,£,iq^y) is 
obtained as 

^ = ^ E E II^Y(z,,.)f - a\M - k) (30) 

where cr^ is expressed in (fT3l) . Hence, a collision status is declared to occur whenever 5 exceeds 
a specified threshold rj. In that case, the BS does not send any response message to the users 
on the considered subchannel as their estimated synch parameters and power levels are regarded 
as not sufficiently reliable. The RSSs that do not find their information repeat the ranging 
process in the next frame using a different ranging opportunity, li 5 < rj, the BS considers 
the users' parameters as accurately estimated and sends a response message for all detected 
codes in the considered subchannel. In the sequel, we refer to the above procedure as the ad- 
hoc collision detector (AHCD). Clearly, r] must be properly designed to achieve a reasonable 
trade-off between the probability of declaring a collision when in fact it is not present (false 
alarm) and the probability of not detecting a collision when in fact it is present (mis-detection). 
Since computing such probabilities by theoretical analysis appears a formidable task, numerical 
simulations can be used in practice for the design of rj. 

In assessing the computational load of AHCD, we \q\. K = K and observe that S{iq^u) is 
available at the BS as an output of AHTE. Hence, computing AY(zg ,^) in (l27l) needs 8KM 
flops for any value of q and u while AM — 1 flops are required to get || AY(ig ,^) ||^. Finally, 
evaluating 6 in (|30l) involves QV real additions. The resulting complexity for the considered 
subchannel is shown in the third line of Tab. I. 
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VI. Simulation results 

A. System parameters 

The investigated OFDMA system is based on the IEEE 802. 16e standard for wireless MANs. 
The DFT size is = 1024 and the sampling period is = 87.5 ns, corresponding to a 
subcarrier distance of \/{NTs) = 11.16 kHz. We assume that A'o = 80 subcarriers are placed 
at both edges of the signal spectrum. The number of modulated subcarriers is thus Nu = 
864 while the uplink bandwidth is approximately = Nu/{NTs) = 9.7 MHz. A total of 
Nji = 144 subcarriers are reserved for ranging. They are divided into -R = 18 subchannels, each 
comprising Q = A tiles uniformly spaced over the signal spectrum at a distance of Nu /Q = 216. 
The number of subcarriers in any tile h V = 2. The remaining Nu — Nji = 720 subcarriers 
are grouped into 15 data subchannels, each composed by 48 subcarriers. A ranging time-slot 
includes M = 4 OFDMA symbols. Hence, the number of ranging opportunities in each frame 
is Ai'totai = -R(M — 1) = 54. The ranging codes are taken from the Fourier basis of order 4, while 
DSS data symbols belong to a QPSK constellation. 

The discrete-time CIRs have maximum order L = 14. Their entries are modeled as independent 
and circularly symmetric Gaussian random variables with zero-mean and an exponential power 
delay profile, i.e., 

^{\hk{i)f} = 4,k ■ exp{-i/Lk), £ = 0, 1, . . . , Lfc - 1 (31) 

where {L^}^^^ are taken from the set {8,9,..., 14} with equal probability while al ^ = (1 — 
e'^)/(l — e"^/^*) so as to have E{||hfe||^} = 1 for all active users. In this way, the powers 
Pk of the received signals are random variables with unit mean. Channels of different users are 
statistically independent of each other. We consider a cell radius of 1.5 km, corresponding to a 
maximum propagation delay of 6'max = 114 sampling periods. Ranging symbols are preceded by 
a CP of length Nq = 128 in order to avoid IBI. The normalized CEOs are uniformly distributed 
over the interval [— ^max, ^max] and vary at each run. Recalling that the estimation range of MEE 
is \£k\ < N / {2MNt), throughout simulations we set ermax < 0.1 for the RSSs. The number 
of CEO trial values in (fT4l) is A^^ = 400, corresponding to an MEE frequency resolution of 
2£max/^e < 5 • 10^''. Eor the DSSs we fix e^ax = 0.02, while the maximum timing error is 
limited to 48 samples. 
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Unless otherwise specified, we consider a static scenario where channel coefficients are 
generated at each simulation run and kept fixed over an entire time-slot. In evaluating the accuracy 
of the synchronization and code detection algorithms, we assume that all RSSs attempt their 
ranging process simultaneously at the first time-slot choosing different ranging opportunities. 
Collision events are simulated only to assess the performance of AHCD. Comparisons are made 
with the ranging scheme proposed by Fu, Li and Minn (FLM) in [8] under a common simulation 
set-up. This includes the same number of ranging subcarriers, ranging subchannels and data 
subchannels, as well as the same transmitted energy from each user terminal. However, as FLM 
can support M different RSS over the same subchannel, the total number of ranging opportunities 
is 72 with FLM and 54 with our scheme. 

In all subsequent simulations, the same number K of RSSs is present in each ranging 
subchannels. This implies that a total of 7r — KR RSSs are simultaneously active in the 
system. Note that letting X = 3 in our ranging scheme corresponds to a fully-loaded system 
where all ranging opportunities are employed. We fix the number of DSSs to 7/5 = 10, although 
extensive simulations indicate that the system performance is only marginally affected by this 
parameter. 

B. Performance evaluation 

1 ) Multiuser code detection: We begin by investigating the performance of MCD in terms of 
mis-detection and false alarm probabilities, say Pmd and Pfa- Fig. 1 illustrates P„td as a function 
of SNR — 1/cr^. The number of active RSSs in each subchannel is either 2 or 3 while the 
maximum normalized CFO is 0.05. The FLM code detector declares the kth ranging code as 
active provided that the quantity 



exceeds a suitable threshold rjpLM which is proportional to the estimated noise power a^. The 
results of Fig. 1 indicate that MCD performs remarkably better than FLM. As expected, the 
system performance deteriorates as K approaches M. The reason is that increasing K reduces 
the noise subspace dimension, thereby degrading the accuracy of MCD. 




Q-l V-1 



(32) 
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Fig. 2 shows Pfa versus SNR for K = \ or: 2 and ^max = 0.05. As is seen, MCD outperforms 
FLM at SNR values greater than of 6 dB. Interestingly, when using FLM the probability of 
false alarm increases with SNR. Such behavior can be explained by observing that the threshold 
f]FLM is proportional to a^, so that it becomes smaller and smaller as the SNR increases. 

2) Frequency estimation: Fig. 3 illustrates the root mean-square-error (RMSE) of the fre- 
quency estimates obtained with MFE vs. SNR for = 2 or 3 and ^max = 0.05. The theoretical 
analysis in (fTT]) is also shown for comparison. As it is seen, the frequency RMSE is approximately 
2 dB worse than its predicted value. The reason is that the result (flTl) is accurate only for 
large values of QV (large data records), while in our simulation set-up QV is limited to 8. 
Again, the system performance deteriorates as K increases. Nevertheless, the accuracy of MFE 
is satisfactory at all SNR values of practical interest. Indeed, an RMSE of 10^^ is obtained even 
with fsT = 3 if SNR > 13 dB. 

The impact of the CEOs on the performance of MFE is assessed in Fig. 4, where the frequency 
RMSE is shown as a function of ^max- The SNR is fixed to 16 dB while K is still 2 or 3. We 
see that the estimation accuracy is only marginally affected by ^max- 

3) Timing recovery: The performance of the timing estimators is assessed by measuring the 
probability P(e) of a timing error event. The latter occurs whenever the estimate O^j"^ gives 
rise to IBI during the data section of the frame. This amounts to saying that the timing error 
A61 = dk dk is larger than zero or smaller than L — Nqd — 1, where Nqd = 48 is the CP 
length during the data transmission period. Note that the mean-shift Ngd/'^ employed in (l23l) 
is also applied to the FLM timing estimator in order to reduce P(e). Fig. 5 illustrates P(e) vs. 
SNR as obtained with AHTE and FLM. The number of active codes in each ranging subchannel 
is K = 2 or 3 while ^max = 0.05. At practical SNR values, we see that AHTE provides much 
better results than FLM. 

In Fig. 6 the timing estimators are compared in terms of their sensitivity to CEOs. For this 
purpose, P(e) is shown as a function of ^max for K = 2 ox 3 and SNR = 16 dB. Again, we see 
that AHTE outperforms FLM, even though the latter is more robust to CEOs. 

4) Power estimation: Fig. 7 illustrates the RMSE of the power estimates as obtained with 
AHPE and the FLM power estimator. In the latter case, the quantity is computed as i\. = 
Zk/{QV) — a'^ I'M with as given in (|32l) . The number of active RSSs in each subchannel 
is either 2 or 3 and Smax = 0.05. The theoretical analysis (t26l) is also shown for comparison. 
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Good agreement between simulation and theory is obtained for K = 2, while a loss of 2 dB is 
observed with K = The reason is that (|26|) has been derived assuming perfect knowledge of 
the frequency offsets, while in practice the accuracy of the CFO estimates degrades with K. At 
low SNR values, both AHPE and FLM provide similar results, but the former takes the lead as 
the SNR grows large. 

5) Impact of channel variations: All previous measurements have been conducted by keeping 
the channel responses fixed during one simulation run. In order to assess the impact of channel 
variations on the system performance, we now consider time-varying channel taps generated 
by passing white Gaussian noise through a third-order low-pass Butterworth filter. The 3-dB 
bandwidth of the filter is taken as a measure of the Doppler bandwidth Bd = vfc/c, where v 
denotes the mobile speed, /c = 2.5 GHz is the carrier frequency and c the speed of light. Figs. 
8 and 9 illustrate the accuracy of MFE and AHTE, respectively, as a function of v with SNR 
= 16 dB. Again, the number of active codes in each ranging subchannel is = 2 or 3 while 
^max = 0.05. The proposed scheme exhibits a remarkable robustness against channel variations 
and can handle a mobile speed of 30 m/s with negligible loss with respect to a static scenario. 

6) Collision detection: Fig. 10 illustrates the performance of AHCD is terms of false alarm 
(Pfa) and mis-detection {Pmd) probabilities as a function of the threshold rj. The false alarm is 
measured in a scenario where two RSSs employing different ranging opportunities are active 
in each subchannel. Vice versa, Pmd is obtained with K = 3 and assuming that two RSSs in 
each subchannel choose the same code. In both situations the SNR is fixed to 16 dB while 
£raax = 0.05. As mentioned in Sect. V, it is important that the probability of mis-detecting a 
collision event is kept as low as possible in order to avoid that RSSs sharing the same ranging 
opportunity adjust their transmission parameters on the basis of incorrect synch information. 
The results shown in Fig. 10 indicate that a Pmd and Pfa of 2 ■ 10"^ can be obtained by setting 
ri = 0.05. 

C. Computational complexity 

It is interesting to compare the proposed ranging scheme with the FLM approach in terms of 
processing requirement in the considered simulation scenario. Assuming that K = 2 RSSs are 
active (on average) in each ranging subchannel, from Table I it follows that 114,000 flops are 
approximately needed by MCD and MFE for each ranging subchannel, while 1,280 flops are 
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required by AHTE and AHPE. The complexity evaluation for FLM is performed in [8], where 
it is shown that code detection involves 1,088 flops while more than 3 Mflops are necessary for 
timing recovery in each subchannel. These results indicate that our scheme allows a significant 
complexity saving with respect to FLM, with a reduction of the overall number of flops by a 
factor 26. 



We have presented a new ranging method for OFDMA systems where uplink signals arriving 
at the base station are plagued by frequency errors in addition to timing misalignments. The 
synchronization parameters of all ranging users are estimated in a decoupled fashion with 
affordable complexity. This is accomplished through a multistage procedure, where user identi- 
fication and CFO estimation is performed first by means of a subspace decomposition approach. 
Frequency estimates are next employed for timing recovery and power level estimation. A simple 
method to detect possible collisions between RSSs employing the same ranging opportunity is 
also investigated. Compared to previous techniques, the proposed scheme exhibits improved 
accuracy with reduced complexity. Computer simulations indicate that the system performance 
is satisfactory even in the presence of frequency errors as large as 10% of the subcarrier spacing. 
The proposed approach can be used to enhance the ranging process of commercial IEEE 802.16- 
based OFDMA systems. 



This Appendix establishes necessary conditions for code identifiability by means of MCD. 
To make the problem analytically tractable, the ranging codes are taken from the Fourier basis 
of order M as expressed in (fT6l) . From (|7]) and (fT6l) we see that matrix C(£:) in ^ has the 
following Vandermonde structure 



VIL Conclusions 



Appendix 
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(33) 
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^j2.[{k-l)/M+e,Nr/N] k=l,2,...,K. (34) 



Recalling that any given code can be univocally identified as long as it is linearly independent 
of all other codes, it follows that C(£) must be full-rank. On the other hand, for a Vandermonde 
matrix the full-rank condition is met if and only if Zk-^ ^ for ki ^ Bearing in mind (|34l) . 
it is easily seen that z^i 7^ ^fc2 is equivalent to putting 



for any integer number i. To proceed further, we reformulate (1351) as 

kfci - ek,\ \Mi + k2 - h\ (36) 

and observe that the right-hand- side cannot be smaller than unity when ki ^ k2. Hence, a 
necessary condition for code identifiability is that MNt l^fci — 1 < 1, which can be ensured 
by setting 

, , N 

kfc < (37) 

' 2MNt 

for k = l,2,...,K. 



References 

[1] "IEEE standard for local and metropolitan area networks: Air interface for fixed and mobile broadband wireless access 
systems amendment 2 : Physical and medium access control layers for combined fixed and mobile operation in licensed 
bands and corrigendum 1," IEEE Std 802.16e-2005 and IEEE Std. 802.16-2004/Cor 1-2005 Std. 2006, Tech. Rep., 2006. 

[2] J. Krinock, M. Singh, M. Paff, A. Lonkar, L. Fung, and C.-C. Lee, "Comments on OFDMA ranging scheme described in 
IEEE 802.16ab-01/0Irl," IEEE 802.16 Broadband Wireless Access Working Group, Tech. Rep., July 2001. 

[3] X. Fu and H. Minn, "Initial uplink synchronization and power control (ranging process) for OFDMA systems," in 
Proceedings of the IEEE Global Communications Conference ( GLOBECOM), Dallas, Texas, USA, Nov. 29 - Dec. 3, 

2004, pp. 3999 - 4003. 

[4] D. H. Lee, "OFDMA uplink ranging for IEEE 802. 16e using modified generalized chirp-like polyphase sequences," in 
Proceedings of the International Conference in Central Asia on Internet (2005), Bishkek, Kyrgyz Republic, Sept. 26 - 29, 

2005, pp. 1 - 5. 

[5] H. A. Mahmoud, H. Arslan, and M. K. Ozdemir, "Initial ranging for WiMAX (802. 16e) OFDMA," in Proceedings of the 
IEEE Military Communications Conference, Washington, D.C., Oct. 23 - 25, 2006, pp. 1-7. 



Mai-ch 17, 2009 



DRAFT 



22 



[6] J. Zeng and H. Minn, "A novel OFDMA ranging methiod exploiting multiuser diversity," in Proceedings of the IEEE Global 

Telecommunications Conference, vol. 1, Washington, D.C., USA, Nov. 26 - 30 2007, pp. 1498 - 1502. 
[7] X. Zhuan, K. Baum, V. Nangia, and M. Cudak, "Ranging enhancement for 802. 16e OFDMA PHY," document IEEE 

802.16e-04/143rl, Tech. Rep., June 2004. 
[8] X. Fu, Y. Li, and H. Minn, "A new ranging method for OFDMA systems," IEEE Transactions on Wireless Communications, 

vol. 6, no. 2, pp. 659 - 669, February 2007. 
[9] M. Wax and T. Kailath, "Detection of signals by information theoretic criteria," IEEE Transactions on Acoustic, Speech 

and Signal Processing, vol. ASSP-33, pp. 387 - 392, April 1985. 
[10] R. Schmidt, "Multiple emitter location and signal parameter estimation," in Proceedings of RADC Spectrum Estimation 

Workshop. Rome Air Development Corp., 1979, pp. 243 - 258. 
[11] M. Morelli, "Timing and frequency synchronization for the uplink of an OFDMA system," IEEE Transactions on 

Communications, vol. 52, no. 2, pp. 296 - 306, Feb. 2004. 
[12] R. Roy, A. Paulraj, and T. Kailath, "ESPRIT - diiection-of-arrival estimation by subspace rotation methods," IEEE 

Transactions on Acoustic, Speech and Signal Processing, vol. 37, no. 7, pp. 984 - 995, July 1989. 
[13] L. Sanguinetti, M. Morelli, and H. V. Poor, "An ESPRIT-based approach for initial ranging in OFDMA systems," To appear 

in IEEE Transactions on Communications, Presented in part at the 9th IEEE Workshop on Signal Processing Advances 

for Wireless Communications (SPAWC), Recife, Pernambuco, Brazil, 2008. 
[14] P. Stoica and R. Moses, Introduction to Spectral Analysis. Englewood Cliffs, NJ: Prentice Hall, 1997. 
[15] P. Stoica and A. Nehorai, "MUSIC, maximum Ukelihood, and Cramer-Rao bound," IEEE Transactions on Acoustics, 

Speech, and Signal Processing, vol. 37, no. 5, pp. 720 - 741, May 1989. 
[16] G. Golub and C. V. Loan, Matrix Computations, Third Edition. Johns Hopkins University Press, 1996. 
[17] L. Sanguinetti, M. Morelli, and H. V. Poor, "A scheme for initial ranging in OFDMA-based networks," in Proceedings of 

the International Conference on Communications (ICC), Beijing, China, May 19 - 23, 2008, pp. 3469 - 3474. 



March 17, 2009 



DRAFT 



23 



TABLE I 

Computational load for each ranging subchannel 



Algorithm 


Required Flops 


MFE & MCD 


AM^iQV + 18M) + 2MN,{4M + 1){M - K) 


AHTE & AHPE 


AK'^M + 2KQV{AM + &K + K'^ + &) 


AHCD 


4:MQV(2K + 1) 




Fig. 1. Mis-detection probability vs. SNR with A" = 2 or 3 and Smax = 0.05. 
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Fig. 2. False alarm probability vs. SNR with K = 1 or 2 and Emax = 0.05. 
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Fig. 3. Frequency RMSE vs. SNR with Ji' = 2 or 3 and Smax = 0.05. 
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Fig. 4. Frequency RMSE vs. £max with A" = 2 or 3 and SNR= 16 dB. 
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Fig. 5. P(e) vs. SNR with if = 2 or 3 and e^ax = 0.05. 
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Fig. 6. P(e) vs. Emax with A' = 2 or 3 and SNR= 16 dB. 
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Fig. 7. RMSE of the power estimates vs. SNR with if = 2 or 3 and Emax = 0.05. 
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Fig. 8. Frequency RMSE vs. v with 7f = 2 or 3 and SNR = 16 dB. 
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Fig. 9. P(e) vs. V with AT = 2 or 3 and SNR = 16 dB. 
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Fig. 10. Performance of AHCD with SNR = 16 dB. 
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